Effects of Preheated Clusters on the CMB Spectrum 
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ABSTRACT 

Mounting evidence from x-ray observations reveals that bound objects should receive some 
form of energy in the past injected from non-gravitaional sources. We report that an instanta- 
neous heating scheme, for which gases in dense regions were subjected to a temperature jump 
of IkeV at z = 2 whereas those in rarified regions remained intact, can produce bound objects 
obeying the observed mass-temperature and luminosity-temperature relations. Such preheating 
lowers the peak Sunyaev-Zeldovich (SZ) power by a factor of 2 and exacerbates the need for 
the normalization of matter fluctuations erg to assume an extreme high value (~ 1.1) for the 
SZ signals to account for the excess anisotropy on 5-arcminute scale detected by the Cosmic 
Background Imager in the cosmic microwave background radiation. 

Subject headings: cosmology: theory - cosmic microwave background - intergalactic medium 



1. Introduction 

In the hierarchical structure formation model 
of cold dark matter, gravitationally bound objects 
ought to obey self-similar scaling relations. That 
is, at any cosmic epoch the mass of a collapse ob- 
ject scales as the third power of its size, and the 
squared velocity dispersion or temperature scales 
as the second power of the size . These scaling re- 
lations follow from a straightforward dimensional 
analysis, given that gravitational collapse occurs 
much faster than Hubble expansion and depends 
on only two dimensional (quasi-)constants: the 
gravitational constant G and the background mat- 
ter density p m Q. Although p m Q varies on the Hub- 
ble time scale, such dependence on long time scale 
affects only the interior profiles of bound objects 
(Navarro et al. 1996; Shapiro & Iliev 2002; Chiueh 
k He 2002). 

However, x-ray observations in the past years 
have consistently indicated that gases in bound 
objects violate the self-similar scaling relations. 
The observed x-ray luminosity-temperature (L — 
T) relation shows that at a given luminosity the 
cluster temperature appears considerably higher 
than the scaling relation (L oc T 2 ), more so 
for low-mass clusters than for high-mass clusters 



(Ponman et al. 1996; Allen & Fabian 1998; Xue 
& Wu 2000). The observed mass-temperature (M- 
T) relation also follows the same trend in violat- 
ing the M oc T 3 / 2 scaling law (Finoguenov et 
al. 2001). More remarkably, recent XMM ob- 
servations of metal lines that directly probe the 
core temperatures of cooling- flow clusters find that 
gases of temperature below some fraction of bulk 
temperature are unexpectedly rare (Peterson et 
al. 2002), suggestive of an entropy floor at the 
cluster core (Ponman, Cannon & Navarro 1999; 
Loewenstein 2000) thus requiring some form of 
non-gravitational energy injection to compensate 
for Bremsstrahlung cooling. 

One candidate for the quested heat sources in 
clusters and groups is the outflow from Active 
Galactic Nucleus (AGN). (Tabor & Binney, 1993; 
David et al. 2001). If this is correct, such a 
mechanism should be more pronounced in the far 
past when quasars were highly active than the 
present. For this and other reasons, investigators 
have tested various schemes for heat injection at 
z > 4 (da Silva et al. 2001; Bialek et al. 2001; 
Muanwong et al. 2002). However, recent com- 
bined analyses of the ASCA data on the entropy 
profiles in galaxy groups and the structure forma- 
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tion model indicate that preheating should occur 
impulsively around z = 2 — 2.5 (Finogucnov et al. 
2002). A similar conclusion has also been reached 
from a different perspective, using the constraint 
on Compton y (< 1.5 x 10~ 5 ) provided by the 
COBE/FIRAS data in testing the SZ effect at the 
hot spots of AGN jets (Yamada & Fujita 2001). 
The present work has been so motivated to con- 
sider impulsive heating around z ~ 2. 

On the other hand, the recent Cosmic Back- 
ground Imager (CBI) experiment discovered that 
the temperature anisotropy of cosmic-microwave- 
background radiation (CMB) exhibits a 3a higher 
level of anisotropy around the harmonic mode 
number 2000 < I < 3500, than anticipated with 
Silk damping (Pearson et al. 2002). This feature 
has been attributed to the large SZ flux (Bond et 
al. 2002). In order to account for such a large SZ 
flux, the value of a& needs to be at least as large 
as unity, though from Figs. 10 and 11 of Bond et 
al.(2002), a slightly larger as fits the CBI data 
even better. 

The SZ flux has been known to depend very 
sensitively on as (Fan & Chiueh 2001; Komatsu 
& Seljak 2002), due to the sensitive dependence 
of cluster number density on the amplitude of 
matter fluctuations at the high end of mass func- 
tion. Hence the SZ effect can serve as a unique 
probe to fix as to high precision. As of today, the 
large-scale galaxy distribution combined with dis- 
tant supernova and CMB data have constrained as 
within a narrow window 0.75 < as < 1 (Bond et 
al. 2002). From the perspective of cluster abun- 
dance explored by the Sloan Digital Sky Survey 
(SDSS) (Bahcall et al. 2002), the extreme value 
as = 1 has, however, placed a considerably low 
ceiling for the matter density fi m = 0.175 ± 0.025. 

If the SZ effect indeed accounts for the CBI 
deep-field excess, preheating can actually exacer- 
bate the awkward situation pertinent to a high 
value of as- Past works have indicated that pre- 
heating generally reduces the SZ flux (Bialcck et 
al, 2001; da Silve et al. 2001). However, these 
works were not able to reproduce both observed 
L — T and M — T relations and hence the reduc- 
tion factor for SZ flux cannot be reliably gauged. 
The present paper reports that impulsive heating 
around z = 2 can reproduce the observed L — T 
and M — T relations and yields a modified SZ 
power spectrum lower than that without preheat- 



ing by a factor 2. It therefore pushes as upward by 
a factor 1.1 for SZ effects to explain the CBI high- 
l result, and fi m downward by 1.2 constrained by 
local cluster abundance. 

2. Simulations, M-T and L-T Relations 

Cosmological hydrodynamical simulations with 
the GADGET code (Springel, Yoshida & White 

2001) were run for both trg = 0.94 and 1 cases in 
a periodic box of comoving size 100/i _1 Mpc on 
our 32-node, dual-cpu Athlon cluster. The simu- 
lations contain 256 3 /2 dark matter particles and 
256 3 /2 gas particles, and the gravitational soft- 
ening length is chosen to be 20h~ 1 kpc. Both 
adiabatic and preheating runs start at z = 100 
with the same initial conditions of n s = 1 and 
(a 8 ,0 TO ,0 6 ,0 A ,/i) = (0.94,0.34,0.05,0.66,0.66) 
and (1,0.3,0.044,0.7,0.64) separately. 

Preheating takes place at z — 2 with a density- 
dependent injection energy per mass, 

Au = u exp[ ^^], (1) 

Pm 

where p m o is the background matter density. Both 
u and P are parameters to be tuned to fit the ob- 
served M — T and L — T relations described be- 
low. This form of energy injection ensures that the 
low-density regions receive no heat, giving rooms 
for some cool gases to survive and the high den- 
sity regions receive a constant energy per mass 
u . The injected entropy has minima at the cores 
of bound objects and maxima at the outskirts 
around p m ~ /3p m o- This heating scheme is sim- 
ilar, though not in detail, to the entropy injec- 
tion of Borgani et al.(2001 & 2002). The best pa- 
rameters are found to be (5 = 1 and m = (500 
km/sec) 2 . Such a value of u corresponds to 1 
keV temperature, and the amount of injected heat 
is equivalent to giving every baryon a 0.5 keV tem- 
perature jump. This amount is roughly consistent 
with those suggested by recent studies on heat in- 
jection from the bound-object interiors (Loewen- 
stcin 2000; Wu, Fabian & Nulsen 2000; Bower et 
al. 2001; Borgani et al. 2001 & 2002; Xue & Wu 

2002) . 

Standard FOF algorithm with a 0.2 linking 
length is employed to locate halos. The mass cen- 
ters may considerably deviate from the peaks of 
spherical over-density, such as in the case of ma- 
jor mergers. In these cases, the over-density so 
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obtained is not useful for testing the L — T and 
M — T relations. They are avoided by selecting 
only those haloes with nonzero M2500, where M n 
represents the mass within a sphere of over-density 
n times that of the critical density. The temper- 
ature T of an ionized plasma is related to the in- 
ternal energy per mass u of the neutral-fluid in 
hydrodynamic simulations by /csT = (2/3)/xm p u, 
where fj, is the mean "molecular" weight of e~, 
p + and He 2+ and equal to 0.588. In addition, we 
define the Bremsstrahlung emission-weighted tem- 
perature (ngT 3 / 2 )/(ngT 1 / 2 ), to be compared with 
the measured temperature of x-ray continuum. 

Figure 1 compares both adiabatic and preheat- 
ing M500-T500 relations from our as = 0.94 and 
1 simulations with the fitting formula given by 
Evrard, Metzler & Navarro (1996) from their adi- 
abatic simulations, and with the fitting formula 
of Finoguenov, Reiprich, & Bohringer (2001) de- 
rived from the observational data re-scaled to z = 
0. Only data with luminosity higher than 10 41 
erg/sec are plotted here. The M 500 — T 50 o relation 
of the adiabatic simulation agrees fairly with that 
of Evrard et a!., but both predict too low a tem- 
perature in bound objects. On the other hand, our 
preheating simulations yield M500 — T500 relations 
of as = 0.94 and 1 agreeing quite well with the 
observations. 

The x-ray luminosity requires somewhat sophis- 
ticated treatments (Mathiesen & Evrard 2001), 
since line emission below 0.5 keV can vastly dom- 
inate the Bremsstrahlung continuum. A constant 
metalicity, 0.3 times the solar value, is adopted 
and the bolometric luminosity reads 

L 50Q = 1.4 x W 27 J T^g^neC^mZiW 

i 

(2) 

in unit of erg/sec, where the volume integration 
is up to i?500j the summation is over different ion 
species, Z is the atomic number, and g(T) the 
Gaunt factor. Figure 2 shows L500 against the 
emission weighted Also shown are the ob- 

served L500 and T500 compiled from Arnaud & 
Evrard (1999), Helsdon & Ponman (2000) and 
Markevitch (1998). The preheating results for 
both as — 0.94 and 1 cases again show excellent 
agreement with x-ray observations. Interestingly, 
Novicki, Sornig & Henry (2002) pointed out that 
the L — T relation is consistent with no evolution 
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Fig. 1. — M500 — T500 relation: Open circles (trian- 
gles) represent adiabatic simulation data for as = 
0.94(1), to be compared with the long-dashed line, 
the simulation results of Evrard et al.(1996). Open 
squares (crosses) are preheating simulation data 
for as = 0.94(1) and are consistent with the solid 
line, the observational result given in Finoguenov 
et al.(2001). All results are for z ~ 0. 

in a low-density universe. Our preheating result 
does display this tendency within < z < 1. This 
behavior also suggests that after initial transients, 
the heated clusters quickly assume certain relaxed 
states since z = 1. 

3. SZ and CMB Power Spectra 

The agreement of our M — T and L — T relations 
with observations confirms the proposed heating 
scheme to be quantitatively viable. We thus pro- 
ceed to construct the 1 deg 2 SZ map by projecting 
the electron pressure through the randomly dis- 
placed and oriented simulation boxes, separated 
by 100/i~ 1 Mpc, along a viewing cone to the red- 
shift direction up to z = 2.5. The average power 
spectra of 40 such SZ maps for each adiabatic and 
preheating simulation of erg = 0.94 and 1 cases are 
shown in Fig. 3. The adiabatic results are consis- 
tent with earlier GADGET results (Springel et al. 
2001), where small-scale structures are more abun- 
dant than those in grid-based hydrodynamic sim- 
ulations (Zhang et al. 2002). Both a s = 0.94 and 
1 runs show that preheating lowers the averaged 
SZ power spectrum by a factor 2 on few-arcminute 
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Fig. 2. — L 50 o — T 500 relation: square, circle, and 
diamond symbols are the observed x-ray luminosi- 
ties and temperatures given by different groups as 
described in the text. The solid, dashed, dotted 
and long-dashed lines are the simulation results 
for (7 g = 0.94 at z =0, 0.5 and 1, and for erg = 1 
at z = 0, respectively. The upper band is for the 
adiabatic simulations, and the lower band for the 
preheating simulations. Both display little evolu- 
tion within < z < 1. 

scales. 

In Fig. 3 we have rescaled the power spectra of 
the cjg = 1 case from original Q^h 2 = 0.018 to 
flbh 2 = 0.022 so as to be compared with those 
analyzed in Bond et al.(2002). The CBI deep- 
field observations scanned a 9 deg 2 sky (Pear- 
son et al. 2002). We hence compute the pre- 
heated SZ power spectrum averaged over 9 one- 
deg 2 SZ maps, which is then added to the (error- 
free) primary CMB spectrum of = 0.04, Q m = 
0.3, f^A = 0.7, n s = 0.975, h = 0.68, parameters 
used in Fig. 11 of Bond et al.(2002), to obtain one 
realization. A total of 40 realizations are con- 
structed. 

The preheated SZ power spectrum of as = 1 
looks very much like the adiabatic power spectrum 
of ag = 0.9, according to the scaling law given by 
Komatsu & Seljak (2001) that the peak SZ power 
approximately scales as a\. The la error bars for 
the as — 1 preheated SZ+CMB power arise from 
the cosmic variance of SZ clusters and are seen to 
be within a factr 1.3 from the mean. Also clear 




Fig. 3. — Comparison of CBI high-! data with the 
averaged preheating and adibatic SZ power spec- 
tra for as = 0.94&1, and with the CMB+(er 8 = 1) 
preheating SZ spectrum. The la error bars of the 
latter arise from the cosmic variance of SZ spectra 
over 40 fields of 9 deg 2 and nearly miss the 2a er- 
ror box of the band-averaged power in CBI deep 
data. 



in Fig. 3 is that the preheated SZ+CMB power 
spectrum nearly misses the 2-a error box of band- 
averaged power in CBI deep data by la. Such 
reduction of SZ power by preheating makes it im- 
plausible for the SZ effect to explain the excess flux 
detected by CBI around 2000 < I < 4000 under 
the framework of concordance model with as = 1. 
Employing the a\ scaling law, we estimate that 
the required as needs to be raised to 1.1 in order 
to be consistent with the CBI high-! excess. 

4. Conclusions 

The agreement of our preheating cosmological 
simulations with the observed M — T and L — T 
relations demonstrates that the impulsive heating 
scheme described by Eq.(l) is a plausible prescrip- 
tion, though it is probably not the unique one and 
its physical origin is unclear. Contrary to our re- 
sults, recent simulations of Borgani et al.(2002) 
adopting entropy injection failed in producing the 
observed M—T relation despite their successful re- 
production of the L — T relation. Various preheat- 
ing schemes so far reported tend to yield reduced 
SZ flux, which in our case can be attributed to that 
gases impulsively blown away from the gravita- 
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tional potential well form an extended halo, which 
can subsequently not fall into, therefore compres- 
sionally heated by, the ever-growing potential well. 
The agreement of our results with the observed 
M—T and L — T relations encourages us to believe 
that our quantitative prediction of the preheated 
SZ flux reduced by a factor 2 is founded on an 
empirically sound base. 

This reduction of SZ flux has a significant impli- 
cation to the interpretation of the CBI deep-field 
results, in that the measured 3cr enhancement of 
CMB high-^ anisotropy is unlikely to be caused by 
the SZ effect of intra-cluster gases in a er 8 = 1 uni- 
verse of concordance model proposed in the orig- 
inal CBI paper (Bond et al. 2002). We estimate 
that as needs to assume a value 1 . 1 for this mech- 
anism to work. The unfortunate consequence of it 
is that Cl m is then pushed to lie about 0.15 ± 0.02, 
constrained by a s ^ 6 = 0.35±0.03 of recent SDSS 
cluster abundance studies (Bahcall et al. 2002). 
Such an £l m can be just too low to be consistent 
with the observed x-ray baryon fraction (< 1/4) 
in massive clusters (Ettori & Fabian 1999; Mohr 
et al. 1999). 
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